257 research outputs found
Symbolic-Numeric Tools for Analytic Combinatorics in Several Variables
Analytic combinatorics studies the asymptotic behaviour of sequences through
the analytic properties of their generating functions. This article provides
effective algorithms required for the study of analytic combinatorics in
several variables, together with their complexity analyses. Given a
multivariate rational function we show how to compute its smooth isolated
critical points, with respect to a polynomial map encoding asymptotic
behaviour, in complexity singly exponential in the degree of its denominator.
We introduce a numerical Kronecker representation for solutions of polynomial
systems with rational coefficients and show that it can be used to decide
several properties (0 coordinate, equal coordinates, sign conditions for real
solutions, and vanishing of a polynomial) in good bit complexity. Among the
critical points, those that are minimal---a property governed by inequalities
on the moduli of the coordinates---typically determine the dominant asymptotics
of the diagonal coefficient sequence. When the Taylor expansion at the origin
has all non-negative coefficients (known as the `combinatorial case') and under
regularity conditions, we utilize this Kronecker representation to determine
probabilistically the minimal critical points in complexity singly exponential
in the degree of the denominator, with good control over the exponent in the
bit complexity estimate. Generically in the combinatorial case, this allows one
to automatically and rigorously determine asymptotics for the diagonal
coefficient sequence. Examples obtained with a preliminary implementation show
the wide applicability of this approach.Comment: As accepted to proceedings of ISSAC 201
Algorithms for Combinatorial Systems: Well-Founded Systems and Newton Iterations
We consider systems of recursively defined combinatorial structures. We give
algorithms checking that these systems are well founded, computing generating
series and providing numerical values. Our framework is an articulation of the
constructible classes of Flajolet and Sedgewick with Joyal's species theory. We
extend the implicit species theorem to structures of size zero. A quadratic
iterative Newton method is shown to solve well-founded systems combinatorially.
From there, truncations of the corresponding generating series are obtained in
quasi-optimal complexity. This iteration transfers to a numerical scheme that
converges unconditionally to the values of the generating series inside their
disk of convergence. These results provide important subroutines in random
generation. Finally, the approach is extended to combinatorial differential
systems.Comment: 61 page
Recursive Combinatorial Structures: Enumeration, Probabilistic Analysis and Random Generation
In a probabilistic context, the main data structures of computer science are viewed as random combinatorial objects.
Analytic Combinatorics, as described in the book by Flajolet and Sedgewick, provides a set of high-level tools for their probabilistic analysis.
Recursive combinatorial definitions lead to generating function equations from which efficient algorithms can be designed for enumeration, random generation and, to some extent, asymptotic analysis. With a focus on random generation, this tutorial first covers the basics of Analytic Combinatorics and then describes the idea of Boltzmann sampling and its realisation.
The tutorial addresses a broad TCS audience and no particular pre-knowledge on analytic combinatorics is expected
On the non-holonomic character of logarithms, powers, and the n-th prime function
We establish that the sequences formed by logarithms and by "fractional"
powers of integers, as well as the sequence of prime numbers, are
non-holonomic, thereby answering three open problems of Gerhold [Electronic
Journal of Combinatorics 11 (2004), R87]. Our proofs depend on basic complex
analysis, namely a conjunction of the Structure Theorem for singularities of
solutions to linear differential equations and of an Abelian theorem. A brief
discussion is offered regarding the scope of singularity-based methods and
several naturally occurring sequences are proved to be non-holonomic.Comment: 13 page
Fast Conversion Algorithms for Orthogonal Polynomials
We discuss efficient conversion algorithms for orthogonal polynomials. We
describe a known conversion algorithm from an arbitrary orthogonal basis to the
monomial basis, and deduce a new algorithm of the same complexity for the
converse operation
Effective Scalar Products for D-finite Symmetric Functions
Many combinatorial generating functions can be expressed as combinations of
symmetric functions, or extracted as sub-series and specializations from such
combinations. Gessel has outlined a large class of symmetric functions for
which the resulting generating functions are D-finite. We extend Gessel's work
by providing algorithms that compute differential equations these generating
functions satisfy in the case they are given as a scalar product of symmetric
functions in Gessel's class. Examples of applications to k-regular graphs and
Young tableaux with repeated entries are given. Asymptotic estimates are a
natural application of our method, which we illustrate on the same model of
Young tableaux. We also derive a seemingly new formula for the Kronecker
product of the sum of Schur functions with itself.Comment: 51 pages, full paper version of FPSAC 02 extended abstract; v2:
corrections from original submission, improved clarity; now formatted for
journal + bibliograph
Efficient Algorithms for Mixed Creative Telescoping
Creative telescoping is a powerful computer algebra paradigm -initiated by
Doron Zeilberger in the 90's- for dealing with definite integrals and sums with
parameters. We address the mixed continuous-discrete case, and focus on the
integration of bivariate hypergeometric-hyperexponential terms. We design a new
creative telescoping algorithm operating on this class of inputs, based on a
Hermite-like reduction procedure. The new algorithm has two nice features: it
is efficient and it delivers, for a suitable representation of the input, a
minimal-order telescoper. Its analysis reveals tight bounds on the sizes of the
telescoper it produces.Comment: To be published in the proceedings of ISSAC'1
On the Complexity of the F5 Gr\"obner basis Algorithm
We study the complexity of Gr\"obner bases computation, in particular in the
generic situation where the variables are in simultaneous Noether position with
respect to the system.
We give a bound on the number of polynomials of degree in a Gr\"obner
basis computed by Faug\`ere's algorithm~(Fau02) in this generic case for
the grevlex ordering (which is also a bound on the number of polynomials for a
reduced Gr\"obner basis, independently of the algorithm used). Next, we analyse
more precisely the structure of the polynomials in the Gr\"obner bases with
signatures that computes and use it to bound the complexity of the
algorithm.
Our estimates show that the version of~ we analyse, which uses only
standard Gaussian elimination techniques, outperforms row reduction of the
Macaulay matrix with the best known algorithms for moderate degrees, and even
for degrees up to the thousands if Strassen's multiplication is used. The
degree being fixed, the factor of improvement grows exponentially with the
number of variables.Comment: 24 page
Fast Computation of Common Left Multiples of Linear Ordinary Differential Operators
We study tight bounds and fast algorithms for LCLMs of several linear
differential operators with polynomial coefficients. We analyze the arithmetic
complexity of existing algorithms for LCLMs, as well as the size of their
outputs. We propose a new algorithm that recasts the LCLM computation in a
linear algebra problem on a polynomial matrix. This algorithm yields sharp
bounds on the coefficient degrees of the LCLM, improving by one order of
magnitude the best bounds obtained using previous algorithms. The complexity of
the new algorithm is almost optimal, in the sense that it nearly matches the
arithmetic size of the output.Comment: The final version will appear in Proceedings of ISSAC 201
- …